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REFINED ALGORITHMS TO COMPUTE SYZYGIES 


BURgiN EROCAL, OLEKSANDR MOTSAIC, FRANK-OLAF SCHREYER, 
AND ANDREAS STEENPASS 


Abstract. Based on Schreyer’s algorithm [S80l IS9ll IBS] . we present 
two refined algorithms for the computation of syzygies. The two main 
ideas of the first algorithm, called |LiftHybrid[ are the following: First, 
we may leave out certain terms of module elements during the compu¬ 
tation which do not contribute to the result. These terms are called 
“lower order terms”, see Definition [472] Second, we do not need to order 
the remaining terms of these module elements during the computation. 
This significantly reduces the number of monomial comparisons for the 
arithmetic operations. For the second algorithm, called [LiftTree] we 
additionally cache some partial results and reuse them at the remaining 
steps. 


1. Introduction 

Computing syzygies, that is, a free resolution 

0 —> F n — > ... —>■ Ei — > Fq — > M — > 0 

of a module M over a polynomial ring R = K[xi,...,x n ], is one of the 
fundamental tasks in constructive module theory, needed for example for 
the computation of the modules Ext^(M, N ) for a further E-module N. An 
algorithm for computing the resolution, starting from a Grobner basis of the 
image of the presentation matrix im(Ei —> Fo), was given by Schreyer [S80| . 
see also [Eis95l Theorem 15.10, Corollary 15.11] or |BS1 Corollary 1.11], 

If M is a graded module over the polynomial ring R. with its standard 
grading, then there exists a minimal free resolution which is uniquely deter¬ 
mined up to isomorphism. While the computation of the Grobner basis of 
the presentation matrix is still feasible in many examples, the computation 
of the minimal free resolution might be out of reach. However, the compu¬ 
tation of a non-minimal resolution is typically much cheaper, and for many 
applications, such as the computation of a single Ext module, good enough. 

In this paper we describe a refined version of Schreyer’s algorithm, which 
utilizes the full strength of Schreyer’s Theorem jEis95i Theorem 15.10]. The 
basic idea, which can be already found in (S91j . is to ignore lower order 
terms in the computation of the generators of the next syzygy module. This 
is possible since these terms will cancel each other anyway. 

Our implementation gives a considerable speed-up, in many cases even if 
we additionally minimize the resolution. The computation of the minimal 
Betti numbers from a non-minimal resolution is typically much faster than 
minimizing the whole resolution. And, for large examples, the computation 
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of a non-minimal resolution using our method is again faster than deriving 
the minimal Betti numbers from it. Note that in many cases, however, even 
a non-minimal resolution suffices to deduce geometric information, see, for 
example, Remark 6.2 

Moreover, our findings suggest that starting, say, from a Grobner basis 
over Q, the computation of a non-minimal resolution using floating point 
numbers as coefficients might be numerically stable. However, this is a 
topic of future work. 

The paper is organized along the following lines. In Section [2j we intro¬ 
duce some basic terminology. The induced monomial ordering, Schreyer’s 
Theorem, and the corresponding algorithm are discussed in Section[3j Based 
on an analysis of this algorithm, we present the two new algorithms in Sec¬ 
tion [4j A detailed example is given in Section [5] In Section [6j we illustrate 
our implementation in a number of examples. One series of examples consists 
of Artinian graded Gorenstein algebras which we regard as an appropriate 
family of examples to test any syzygy algorithm since we can vary the num¬ 
ber of variables, the degree, and the sparseness. Further series are nodal 
canonical and nodal Prym canonical curves. For the examples, we always 
work over a finite ground field. 


2. Preliminaries 


Throughout this article, let K be a field, and let R := K[x i,..., x n \ be the 
polynomial ring in n variables over K. We denote the monoid of monomials 
in xi,...,x n by Mon(xi,.. .,x n ). 

We briefly recall some terminology for dealing with R-module syzygies 
and their computation. 

Definition 2.1. Let F := R r be the free 17,-module of rank r, and let 
ei,..., e r be the canonical basis of F. 

(1) A monomial in F is the product of an element in Mon(xi,... ,x n ) 
with a basis element e*. The set of monomials in F is denoted by 
Mon(F’). 

(2) Accordingly, a term in F is the product of a monomial in F with a 
scalar in K. 

(3) A monomial miej divides a monomial if i = j and mi divides 
m 2 ; in this case, the quotient m 2 ^j /m\ei is defined as m 2 /mi G 
Mon(xi,..., x n ). We also say that mi G Mon(xi,..., x n ) divides 
rri^dj if mi divides m 2 , and in this case m^e^/mi is defined to be 
(m2/mi)ej G Mon(F). 

(4) The least common multiple of two monomials miej,m 2 ej G F is 



(5) A monomial ordering on F is a total ordering >- on Mon(T) such 
that if miej and m^e^ are monomials in F, and m is a monomial in 
R, then 


miej >~ m, 2 ej => (m • mi)ej >~ (m • ni 2 )ej . 
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In this article, we require in addition that 

m\Qi >~ m, 2 &i <=> miej y rn^ej for all i,j. 

(6) Let ^ be a monomial ordering on F, let / E F \ {0} be an element 
of F. and let / = crnet + f* be the unique decomposition of / with 
c E K \ {0}, me* E Mon(T), and mej > m*ej for any non-zero 
term c*m*ej of /*. We define the leading monomial, the leading 
coefficient, the leading term, and the tail of / as 

LM(/) := met , 

LC(/) := c, 

LT(/) := cmei, 
tail(/) := / - LT(/), 

respectively. 

(7) For any subset S C F, we call 

L(S) := (LM(/) | feS\{0}) R cF 
the leading module of S. 


Remark 2.2. Let A be a monomial ordering on the free R -module F := R r 
as defined above. Then there is a unique monomial ordering > on R which 
is compatible with >- in the obvious way, and we say that >- is global if > is 
global. In this article, all monomial orderings are supposed to be global. 


Definition 2.3. Let N be an 77-module, let G := {/i,..., f r } C N be a 
finite subset of N, and let F := R r be the free R- module of rank r as above. 
Consider the homomorphism 

ipG '■ F N, 
e* H- fi. 

A syzygy of G = {/i,..., f r } is an element of ker ifc- We call ker ipc the 
(first) syzygy module of G, written 

Syz(G) := ker ifc- 


Definition 2.4. Let M be an 72-module. A free resolution of M is an exact 
sequence 

F: ... —> F i+1 ^4 Ft A F 4 _i —> ... —> Fi A F 0 —> M —> 0 
with free 72-modules Fj, i G N. 


Remark 2.5. Let the notation be as in Definitions 12.31 and !2. 41 In this arti¬ 
cle, we only consider the case where N is a free module over the polynomial 
ring 72 and where we wish to construct a free resolution of M = N/(G)r. 
For this, with notation as in Definition 2.3, set F$ := N, F\ := F, and 


4>i '■= f’G- Now, starting with G\ := G, let G{+\ be a finite set of generators 
for Syz(G'j) and, inductively, define <f>i to be the map fjc, for i E N\{0}. We 
then have Syz (Gf) = ker (fi, that is, F is obtained by repeatingly computing 
the syzygies of finite subsets of free 72-modules. 
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Definition 2.6. Let Fq := R s be the free R-module of rank s, let > be a 
monomial ordering on Fq, and let G := {/i,..., f r } C Fq \ {0} be a set of 
non-zero vectors in Fq. 

(1) We define rriji as 


rriji 


LCM(LM(/j), LM(/j)) 
LT(/0 


€R. 


(2) For i,j G r - }, we define the S-vector of j) and fj as 


S(fi, fj) := rrijifi ~ mjfj € (G) R C F 0 . 


(3) For g G Fq, we call an expression 


9 = 9ifi + • • • + 9rfr + h 


with gt G R and h G Fj a standard representation for g with remain¬ 
der h (and w.r.t. G and >) if the following conditions are satisfied: 

(a) LM(< 7 ) > LM (gifi) for alH = 1,..., r whenever both g and gift 
are non-zero. 

(b) If h is non-zero, then LT (h) is not divisible by any LT(/j). 


Remark 2.7. Standard representations can be computed by multivariate 
division with remainder. With notation as above, let now G be a Grobner 
basis, and let g be an element of ( G)r. In this case, the remainder h is zero 
by Buchberger’s criterion for Grobner bases. For S-vectors of elements of G , 
each standard representation 

S(fi, fj) = rriji fi - rriijfj = g[ lj) fi + ... + f r 
yields an element m^e* — niijOj — (g^e i + ... + gr^e r ) G Syz(G). 

This gives one possibility to compute syzygies which we will now discuss 
in detail. 


3. Schreyer’s Syzygy Algorithm 

3.1. The Induced Ordering. 

Definition 3.1. Given a monomial ordering > on Fq := R s and a set of 
non-zero vectors G := {/i, ..., f r } C Fo\{0}, we define the induced ordering 
on F\ ■= R r (w.r.t. > and G) as the monomial ordering >- given by 

m\ei >- m 20 j LT(rrti /,) > LT(rr^fj) 

or (LT(mi/j) = LT (m 2 fj) and i > j) 

for all monomials mi, m 2 G Mon(an ..., x n ), and for all basis elements 
ei,ej G F 1 . 

This definition implies that both > and >- yield the same ordering on R 
if restricted to one component. 

Monomial comparisons w.r.t. induced orderings are computationally ex¬ 
pensive and should therefore be avoided in practice. This holds in particular 
in the case of chains (>~i)i= of orderings with Fj+i induced by Fi which 
appear in the computation of free resolutions. 
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3.2. Schreyer’s Theorem. 

Theorem 3.2 ( |BSl Corollary 1.11]). 

Let G = {/i,-..,/r} C Fq := R s be a Grobner basis w.r.t. a monomial 
ordering > on F$. For each pair (/$, fj ) with i,j E {1,... ,r}, let 

S(fi, fj) = rrijifi - mijfj = g^ fi + ■ ■ ■ + 9^ fr 

be a standard representation of the corresponding S-vector. Then the rela¬ 
tions 

rrijifi ~ mijZj - (g^e i + ... + g^e r ) E F, := IF 

form a Grobner basis of Syz(G) w.r.t. the monomial ordering on F\ in¬ 
duced by > and G. In particular, these relations generate the syzygy module 
Syz(G). 

Based on this theorem, there is an obvious algorithm for the computa¬ 
tion of syzygy modules: Given a Grobner basis G as above, it suffices to 
compute standard representations for all S-vectors S(fi, fj) by division with 
remainder. 

Of course, one can do much better. Since S(fi,fj ) = —S(fj,fi), it is 
sufficient to consider those pairs ( fi,fj ) with j < i. It is well-known that 
even more pairs can be left out using the following notation (cf. [S3]): 

Notation 3.3. Let Fo := R s be the free i?-module of rank s, and let G := 
{/i,..., f r } C Fq \ {0} be a set of non-zero vectors in Fo. For i = 2, ..., r, 
we define the monomial ideal Mi as 

Mi := (LT(fi ),... ,LT(/ 4 _!)) : (LT(/j)) C R. 


Remark 3.4. Recall that if N\ and N 2 are submodules of an F-module M, 
then the module quotient N\ : N 2 is defined to be the ideal 


N\ : N 2 := {a E R \ an E N\ for all n E N 2 } C R . 


In particular, in the situation of Notation 3.3, we have (mie'f) : (m^e'-) = 0 
for any two monomials m\ , m 2 E R and any two basis elements e [, e'- of Fo 
with i j. 


Proposition 3.5 (IBS, Theorem 1.5]). Let G = {/i,...,/ r } C Fo be as 
in Theorem \3.H\ For each i = 2,... ,r, and for each minimal generator x a 
of the monomial ideal Mi C R, let j = j(i,a) < i be an index such that 
niji divides x a . Then it is sufficient in Theorem 3.2 to consider only the 
corresponding pairs (fi, fj). 

Taking this proposition into account, we get Algorithm [l] below. 


3.3. Schreyer Frame. The leading module of the syzygy module will serve 
as a starting point for the algorithms which we propose in Section [4] Its 
computation is based on the following observation. 


Remark 3.6. With notation as in Theorem 


3.2 


9l j) fl + • • • + 9r j) fr 


IS 


a standard representation of the S-vector S(fi,fj), and therefore we have 
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Algorithm 1 SyzSchreyer 

Input: A Grobner basis G = {/i,..., f r } C Fq := R s w.r.t. some monomial 
ordering > 

Output: A Grobner basis of Syz(G) C F\ := R r w.r.t. the monomial or¬ 
dering induced by > and G 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


S :=0 

for i = 2 , ..., r do 

for each minimal generator x a of the monomial ideal Mi do 
choose an index j < i such that niji divides x a 
h ■= S(fi, fj) = nrijifi - rriijfj G F 0 
s := viji^i — niijej G F\ 

while h 7 ^ 0 do 

choose an index A such that LT(/\) divides LT(/t) 

LT th) 

s s LT ( /A) e a 
5 := SU{s} 

return S 


LM {rrijifi) = LM (rriijfj) > LM (g^ fk) for all k = 1,. .. ,r with / 0, 
cf. Definition 2.6[ 3 ). For i > j, this implies 

mjiei >- mijej y LM {g^ j) ^je k , 

where >- is the monomial ordering on F\ = R r induced by > and G. There¬ 
fore the leading syzygy module of G w.r.t. >- is 


L^(Syz(G))= © 


Mi 


z=2,...,r 

Thus, for a given Grobner basis G, the leading module of Syz(G) w.r.t. 
the induced ordering can be easily computed by throwing away superfluous 
elements, see Algorithm [2} 

In Algorithm [2j only the leading terms of the Grobner basis G contribute 
to the computation of the set S (via the term rriji G R, cf. Definition PPh 
For a free resolution as constructed in Remark |2.5[ we can thus, starting with 
the leading terms of G, inductively compute sets of generators for all leading 
syzygy modules. The sequence of these sets of leading syzygy terms is called 
a Schreyer frame by La Scala and Stillman in [L5]. 

It is worth noting that the algorithm to compute a minimal free reso¬ 
lution by La Scala and Stillman is compatible with our algorithms for the 
computation of syzygies in the sense that both approaches are based on the 
Schreyer frame and can thus be combined. 

Remark 3.7. In the computation of a free resolution, reordering the syzy¬ 
gies after each step may yield smaller generators for higher syzygy modules. 
With notation as in Remark 2.5, we expect that reordering Gi w.r.t. the 
negative degree reverse lexicographical ordering on Fi-\ before computing 
Gi +1 is generally the best choice. 
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Algorithm 2 LeadSyz 

Input: A Grobner basis G = {/i,..., f r } C Fq := R s w.r.t. some monomial 
ordering > on Fq 

Output: A minimal set £ of generators for the leading syzygy module 
L^(Syz(G)) of G w.r.t. the monomial ordering >- on F\ := R r induced 
by > and G 

1 : C := 0 

2: for 1 < j < i < r do 
3: t := rrijiei £ F\ 

4: for s £ C do 

5: if s | t then 

6: t := 0 

7: break 

8: else if t \ s then 

9: C := C \ {s} 

10: if t 7 ^ 0 then 

11: £:=£U{t} 

12: return C 


4. New Algorithms 


Throughout this section, let G := {/i,..., f r } C Fo : = F s be a Grobner 
basis w.r.t. some monomial ordering > and let >- be the monomial ordering 
on F\ := R r induced by > and G. Furthermore, let C be the minimal gen¬ 
erating set of the monomial submodule L^(Syz(G)) C F\. We simply write 
for the map i/’g : F\ —> Fq dehned by i/jq (e ? ;) := fi as in Definition 2.3 


By Remark 3.6, there is a one-to-one correspondence between the mini¬ 
mal generators of the monomial ideals M t and the elements of C. Instead of 
processing S-pairs, we can therefore directly start with the minimal gener¬ 
ating set of leading syzygy terms. This is equivalent to applying the chain 
criterion for syzygies to the set of all S-pairs, cf. IGPi. Lemma 2.5.10]. 

The algorithmic idea is that each leading syzygy term s £ C gives rise 
to a pair of indices (i, j) with s = rrijiei , which, through a standard repre¬ 
sentation of the corresponding S-vector S(fi, fj), gives rise to a syzygy s of 
G with LT^(s) = s. Note that both the pair of indices and the standard 
representation obtained thereof are in general not unique. 

This motivates the following definition. 


Definition 4.1. Let s £ L^(Syz(G)) C F\ be a leading syzygy term. We 
call s £ F\ a lifting of s w.r.t. G and >- if the following conditions hold: 

(1) LT^(s) = s, and 

(2) s e Syz(G'). 

If we know how to compute such a lifting, then we can use Algorithm [3] 
to obtain a generating set S of the syzygy module. Since L^(5) is equal 
to L^(Syz(G)), this set is even a Grobner basis of Syz(G) w.r.t. >-. From 
the computational point of view, Algorithm [l] can be regarded as the special 
case of Algorithm [3] where the liftings are computed by the usual reduction. 
This can be reformulated as in Algorithm [4| 
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Algorithm 3 SyzLift 

Input: A Grobner basis G C Fo w.r.t. > and an algorithm Lift to compute, 
for a leading syzygy term s £ L^(Syz(G)), a lifting w.r.t. G and A 
(Lift can be, for example, any of the three algorithms LiftReduce, 
LiftHybrid, or LiftTree below.) 

Output: A Grobner basis of Syz(G) C F\ w.r.t. A 

1: C := LeadSyz(G) 

2: S := 0 
3: for s £ C do 
4: s := Lift(s) 

5: S:=SU{s} 

6: return S 


Let us now discuss algorithms for lifting leading syzygy terms in detail. 
LiftReduce (Algorithm |4| computes a lifting of a given leading syzygy 
term s £ L^(Syz(G)) via multivariate division of the polynomial g := ip(s) £ 
(G) C .Fo by the elements of G. This is computationally the same as the 
division of h w.r.t. G in the while-loop of SyzSchreyer (Algorithm [TJ . At 
each step, the leading term of g is reduced, and this process finally reaches 
g = 0 since G is a Grobner basis. 


Algorithm 4 LiftReduce 

Input: A Grobner basis G = {/i,...,/ r } C Fq w.r.t. > and a leading 
syzygy term s £ L>_(Syz(G)) C F x 
Output: A lifting s £ Syz(G) C F\ of s w.r.t. G and >- 

1: g := ijj(s) 

2: s := s 

3: while g ^ 0 do 

4: t := LT( 5 ) 

5: choose a term me* £ F\ with mLT(/j) = t and s A me, 

6: g ■= g- mfi 
7: s := s — mei 

8 : return s 


Let g \,..., gk £ Fo be the sequence of values which g takes when the al¬ 
gorithm LiftReduce is applied to a leading syzygy term s £ L^(Syz(G)). 
Since we have gk = 0, every single term occurring in this sequence is even¬ 
tually cancelled at one of the reduction steps in line [6j but only the pro¬ 
cessing of the leading terms LT(gi),..., LT(^fc) contributes to the syzygy 
s £ Syz(G). In particular, those terms which are not divisible by one of 
the leading monomials LM(/j), i = 1,..., r, do not contribute to s and can 
therefore be left out. We use the following terminology to refer to these 
terms. 

Definition 4.2. Let S C Fo be a set of vectors and let t £ Fo be a term. 
Then t is called a lower order term w.r.t. S if 

LM(/) \t for all / £ S \ {0} . 
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For an element g G Fq, we define LOTUS') to be the sum of those terms 
occuring in g which are of lower order w.r.t. S. 

Furthermore, instead of reducing the leading term of g at a given step, 
we may choose any term of g which is not of lower order. Taking the above 
observations into account, we get the algorithm LiftHybrid (Algorithm [5]). 
Note that the lower order terms which are left out at the intermediate steps 
sum up to zero. 

Algorithm 5 LiftHybrid 

Input: A Grobner basis G = {fi, ■ ■ ■, f r } C Fq w.r.t. > and a leading 
syzygy term s G L>_(Syz(G)) C F\ 

Output: A lifting s G Syz(G) C F\ of s w.r.t. G and 

1: g:=^(s)-LOT^(s)\G) 

2: s := s 

3: while g ^ 0 do 

4: choose a term t of g 

5: choose a term me* G F\ with rn LT(/j) = t and s >- mei 

6: g := g - (mfi-LOT(mfi\G)) 

7: s := s — mei 

8: return s 


We can even go further and consider the set T of terms in g rather than 
the polynomial g itself. In other words, we do not need to sort the terms in g 
and we do not need to carry out the cancellations of terms which may occur 
in line[6]of LiftHybrid. Then each term in T can be reduced independently 
as in LiftTree (Algorithm [6]) . This yields a tree structure by the recursive 
calls of LiftSubtree (Algorithm [6a| for each term in T. 

The algorithm applied at the root node of this tree, LiftTree, slightly 
differs from the algorithm applied at the other nodes, LiftSubtree. In 
LiftTree, the leading term of if(s) is included in T, whereas at the other 
nodes, this term has been cancelled by the reduction in the previous step 
and is therefore left out in LiftSubtree. Because of this difference, we 
need the following definition to give a proper description of the output of 
LiftSubtree. 

Definition 4.3. Let s G F\ be a term. We call s G F\ a subtree lifting of s 
w.r.t. G and if the following conditions hold: 

(1) LT^ (s) = s, and 

(2) all terms in tail(^>(s)) G Fq are lower order terms w.r.t. G and >-. 

LiftTree terminates when T = 0 is reached in every branch of the tree. 
One can easily check that this algorithm returns indeed a lifting of the input 
by comparing it to LiftHybrid. 

Remark 4.4. Let s G L^(Syz(G)) C F\ be a leading syzygy term. If s is a 
lifting of s w.r.t. G and then s is a subtree lifting of s, but the converse 
statement is not true in general. 

For any term s' G F \, a proper lifting of s' (w.r.t. G and >~) exists if and 
only if s' is leading syzygy term, that is, an element of L^(Syz(G)). Hence 
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Algorithm 6 LiftTree 

Input: A Grobner basis G = {/i,...,/ r } C Fq w.r.t. > and a leading 
syzygy term s £ L>_(Syz(G)) C F\ 

Output: A lifting s £ Syz(G) C F\ of s w.r.t. G and >- 

1: g := ip(s) 

2: T := set of terms in (g — LOT^G)) 

3: s := s 

4: for all t £ T do 

5: choose a term me* € Fj with rn LT(/j) = t and s y mei 

6: s := s - LlFTSUBTREE(mej) 

7: return s 


Algorithm 6a LiftSubtree 

Input: A Grobner basis G = {/i,. .., f r } C To w.r.t. > and a term s £ F± 
Output: A subtree lifting s £ F\ of s w.r.t. G and >~ 

1: g := ip(s) - LT(i/>(s)) 

2: T := set of terms in (g — LOT(g r |G ! )) 

3: s := s 

4: for all t £ T do 

5: choose a term me* £ F\ with rn LT(/j) = t 

6: s := s — LlFTSUBTREE(mej) 

7: return s 


we cannot expect to find proper liftings of the terms mei £ F± to which 
LiftSubtree is applied in Algorithms [6] and |6a[ 

Remark 4.5. The condition s me, in line [5] of Algorithm [6] is always 
satisfied at the analogous step in Algorithm [6a] and thus does not need to 
be checked there. 


Remark 4.6. Since the only differences between the algorithms LiftTree 
and LiftSubtree are the assignment of g in line 1 and, as explained in 
Remark 4.5 above, the condition s y met in line 5, we could have merged 
them both into one algorithm. This can be done, for example, by using a 
boolean variable which is set to true in the case of LiftTree, corresponding 
to the root node of the resulting tree, and which is set to false for LiftSub¬ 
tree, representing the inner nodes and the leaves. However, we think that 
separating the two algorithms may help to understand the mathematical 
properties of the output of LiftSubtree as described in Definition 4.3 in 
contrast to the output of LiftTree, see Definition 4.1 


Remark 4.7. Let mie^,... ,mkei k £ F\ be the sequence of terms chosen 
in line [5] of LiftReduce when this algorithm is applied to a leading syzygy 
term s £ L>_(Syz(G)). Then we have 


S y mie h y ... y m k e ik . 

However, the terms mei £ F± chosen in LiftHybrid, LiftTree, and Lift¬ 
Subtree satisfy s y mei , but they are not necessarily ordered. 
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LiftTree has two main advantages in comparison to LiftHybrid. First, 
no reductions as in line [6] of LiftHybrid occur. Second, the results of 
LiftSubtree can be cached and reused. We will see an example for this in 
the next section. 


Remark 4.8. It is worth mentioning that the proposed algorithm can be 
easily parallelized. First of all, it is inherently parallel in several ways: 
Whenever a single leading syzygy term is lifted to a syzygy by LiftTree, 
the different branches of the resulting tree, which correspond to the recursive 
calls of LiftSubtree, are independent of each other and can thus be treated 
in parallel. (Note that this approach should be implemented in such a way 
that it works well with the caching of partial results.) Likewise, in the 
computation of a syzygy module via Algorithm [3j the leading syzygy terms 
can be treated in parallel. In view of a whole resolution, we may start to 
lift the leading syzygy terms in the lower syzygy modules while we are still 
computing the Schreyer frame, see Section 3.3, for the higher ones. Note 
that, however, the time to compute the Schreyer frame is almost negligible 
in most cases. 

Over the rational numbers, one could also apply modular methods to our 
algorithm, that is, one could do the computation modulo several primes in 
parallel and then lift the results back to characteristic zero. Both approaches 
are, however, subject to future research. 


5. Example 

In this section, we give an example in order to illustrate the differences 
between the three approaches which we presented in the previous section. 
The example has been chosen in such a way that it shows the benefits, 
but also possible drawbacks of the new methods. However, note that a 
considerable speed-up can only be expected for large examples. 

Throughout this section, let Fq := R := Q[w, x, y, z] be endowed with 
the lexicographical ordering, denoted by >. We compute the first syzygy 
module of G := (/i, / 2 , / 3 ) C R with 

fi := wx + wz + x 2 + 2 xz — z 2 , 

/2 := wy — wz — xz — yz — 2 z 2 , 

h ■= xy + z 2 . 

Note that G is a Grobner basis w.r.t. >. Let >»- be the Schreyer ordering 

on Ei := R 3 induced by > and G. We can use Algorithm [ 2 ] to check that a 

minimal generating set of the leading syzygy module of G w.r.t. >- is given 
by 

£ = {x ■ e- 2 , w • 63 } C F\ . 

Our goal is to extend these leading syzygy terms to generators of the 
syzygy module. Let us first consider the usual LiftReduce approach (Al¬ 
gorithm [4]). Flow charts of LiftReduce applied to the two leading syzygy 
terms above are shown in Figure [lj 

They start with the input term on the syzygy level and its image g under 
■0 : F\ —> Fq, e* i-A fi, on the level of Fq. At each step, the leading term of 
g is reduced w.r.t. G while si and S 2 keep track of these reductions. Both 
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Figure 1. LiftReduce/LiftHybrid applied to 
x ■ e 2 and w ■ e 3 



charts have the shape of a chain because every step depends on the previous 
one. We could choose a different reduction at the first step of the diagram 
on the right hand side, but there is no other choice at the other steps. The 
process ends when g = 0 is reached, and we finally get the syzygies 

ai = (—y + z) ■ ei + (x + z) ■ e 2 + (x + 3z) • e 3 G Syz(G) and 
S 2 = ~y ■ ei + z • e 2 + (w + x + 2z) ■ e 3 G Syz(G') 

as liftings of the leading syzygy terms x ■ e 2 and w ■ e^, respectively. 

The main innovation of the algorithm LiftHybrid (Algorithm [ 5 ]) is to 
leave out the lower order terms in g. This in turn allows us to choose, at each 
step, any of the remaining terms for reduction, in contrast to LiftReduce. 
Hence the terms in g do not have to be ordered at all. If we always choose, 
however, to reduce the leading term as in LiftReduce, then the flow charts 
of LiftHybrid applied to x-e 2 and u;-e 3 , respectively, can be obtained from 
those for LiftReduce by leaving out the underlined lower order terms, see 
Figure [lj 

In LiftTree (Algorithm [6]), the polynomial g is replaced by a set of terms 
denoted by T and each term is treated independently. The corresponding 
flow charts in Figure [2] and Figure [3] thus have a tree structure where each 
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Figure 2. LiftTree applied to x ■ e 2 



Figure 3. LiftTree applied to w ■ 



node represents one of the recursive calls of LiftTree and LiftSubtree. 
In Figure [3j the result of LiftTree(® • e 2 ) can be read off as the sum of all 
the terms si. Similarly, LiftTree (w- e^) yields — y-ei+z-e2+(w+x+2z)-e3. 

Again, the underlined lower order terms are left out. The process ends 
when T = 0 is reached in every branch. It is worth noting that although each 
step resembles a reduction step, no reductions as in the first two approaches 
occur. The main advantage of the LiftTree approach is that intermediate 
results can be cached and reused. In Figure [3j the term wxy occurs as an 
element of T, but the whole subtree which corresponds to this element has 
already been computed when LiftTree was applied to x ■ e 2 in Figure [2] 
and we can therefore just plug in the cached result. 

A possible drawback of this method can be observed in Figure [2j Two 
steps are necessary to compute the term (3 z ■ e 3 ) in the result whereas 
LiftReduce and LiftHybrid need only one step for this. On the other 
hand, we could also cache the result of LiftSubtree(z • 63 ) and reuse it for 
the computation of LiftSubtree(2z • 63 ), of course. 

As mentioned above, the leading terms of the two computed syzygies 

si = (—y + z) ■ ei + (x + z) ■ Q 2 + (x + 3 z) ■ e 3 G Syz(G) and 
S 2 = -y ■ ei + z ■ e 2 + (w + x + 2z) ■ e 3 G Syz(G) 
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w.r.t. the Schreyer ordering are x ■ e 2 and w ■ e 3 , respectively. Since these 
two terms belong to different module components, the S-vector S(s\,S 2 ) of 


the two syzygies is 0. This implies Syz(Syz(G)) = 0 by Theorem 3.2 applied 
to Syz(G). We therefore get 


T 


0 


R 


2 ^ 2 -> R 3 — 


R 


R/(G) 


R 


0 


as a free resolution of R/(G)r , where the maps are given by := fi, for 

i = 1,2,3, and 02 (?]) '■= Sj, for j = 1,2. Since G and {si,S 2 } are reduced 
Grobner bases w.r.t. the monomial orderings > and respectively, J- is 
even a minimal free resolution. 


6 . Timings and Statistics 


In this section, we illustrate, by a number of examples, the speed-up 
achieved by LiftTree in comparison to other algorithms, and we also give 


detailed statistics on these computations. In Subsection 6.1 we consider two 


series of Artinian graded Gorenstein rings while Subsection 6.2 is devoted 


to randomly constructed canonical and Prym canonical nodal curves. 

For each example, we compare the timings to compute a free resolution 
by the Macaulay2 fM2j command res, by the Singular [ DGPSl command 
lres(), and by the implementation of LiftTree in the Singular library 
schreyer.lib |M14| . Note that the first two yield minimal free resolution 
whereas the resolutions computed by LiftTree are in general non-minimal. 
We therefore additionally present the timings to compute the minimal Betti 
numbers via the Singular command bettiO and to minimize the whole 
resolution via the Singular command minresO, starting from the out¬ 
put of LiftTree in both cases. In separate tables, we also list, for each 
example, the number of terms in the resolution computed by LiftTree 
(excluding the first map which is given by the input ideal) as well as the 
number of multiplications, additions, and cancellations (that is, additions 
to zero) of coefficients in the ground held which were performed during this 
computation. To measure the sparseness of the resolution, the quotient of 
the number of terms divided by the number of (matrix) entries Q, 


sparse 


IS 


also given. Finally, we present the (minimal or non-minimal) Betti tables 
for selected examples. 

In each case, we used the degree reverse lexicographic monomial ordering 
(dp in Singular) and we computed a reduced Grobner basis of the input 
ideal w.r.t. this ordering beforehand. The timings were computed on an Intel 
Core i7-860 machine with 16 GB RAM and 4 physical (8 virtual) cores, each 
with 2.8 GHz, running Fedora 20 (Linux kernel version 3.17.4). A dash (—) 
indicates that the computation did not finish within 24 hours. 


6.1. Artinian graded Gorenstein rings. Our first series of examples con¬ 
sists of Artinian graded Gorenstein rings (AGR) in n + 1 variables of socle 
degree d (see Tables [l] to [9]). We regard this family of zero-dimensional ideals 
as a family where we can easily alter the number of variables, the degree, and 
the sparseness. Of course it would be nice to produce higher dimensional 
random examples, say, random ideals of codimension c < n generated by s 
forms of degree d\,... ,d s with s ^ c. However, since the Hilbert schemes 
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of such examples are not unirational in most cases, we do not know how to 
construct such examples. 

Artinian graded Gorenstein rings in n + 1 variables of socle degree d can 
be obtained from a homogeneous form / on P n of degree d via apolarity, 
see, for example, |RS] . If we write / = + ... + i d s as a sum of s powers of 

linear forms, then the minimal such s is called the Waring rank of /. The 
shape of the minimal resolution and even the length of the Gorenstein ring 
depends on s, and we invite the reader to prove in our examples that for 
small s , the Waring decomposition is unique up to scalar, a topic already 
studied for binary forms by Sylvester and recently picked up again in [ 00 ]. 

Below we report the findings for random examples over the finite field 
with 10,007 elements. It requires 24.12 seconds to compute the minimal 
resolution for s = 48, while it just takes 3.12 seconds to get the minimal 
Betti numbers. 

The total number of additions in most of these examples is smaller than 
the total number of terms, and the number of cancellations (additions to 
zero) is even much smaller. Furthermore, the matrices in the non-minimal 
resolution have basically monomial entries as Figure [4] indicates. This sug¬ 
gests that the algorithm LiftTree to compute a non-minimal resolution is 
numerically stable. 

Table 1. Timings for the AGR examples with d = 5 and 
n = 6 (in sec.) 


s 

Macaulay2 

Singular 

res 

Ires() 

LiftTree 

betti() 

minres() 

12 

5.28 

6.35 

0.06 

0.02 

0.96 

18 

10.65 

3.63 

0.27 

0.06 

4.53 

24 

26.89 

7.59 

1.06 

0.19 

33.31 

30 

16.77 

5.06 

0.94 

1.61 

39.26 

36 

44.45 

28.16 

0.73 

2.09 

11.74 

42 

87.65 

48.85 

0.73 

2.38 

23.39 

48 

87.61 

48.69 

0.73 

2.39 

23.39 


Table 2. Statistics for the non-minimal resolutions of the 
AGR examples with d = 5 and n = 6 


s 

sec. 

# Terms 

#Mult. 

#Add. 

#Canc. 

Qsparse 

12 

0.06 

34,963 

74,190 

38,312 

2,163 

0.155 

18 

0.27 

123,144 

700,889 

573,143 

7,267 

0.224 

24 

1.06 

316,492 

5,961,627 

5,638,864 

16,813 

0.276 

30 

0.94 

319,580 

627,508 

315,310 

2,496 

0.177 

36 

0.73 

294,730 

447,245 

162,996 

331 

0.195 

42 

0.73 

294,762 

447,249 

163,002 

324 

0.195 

48 

0.73 

294,746 

447,260 

162,992 

334 

0.195 
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Table 3. Minimal Betti table for the AGR example with 
n = 6, d = 5, and s = 18 



0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

- 

- 

- 

- 

- 

- 

- 

1 

- 

10 

4 

- 

- 

- 

- 

- 

2 

- 

- 

60 

136 

130 

60 

11 

- 

3 

- 

11 

60 

130 

136 

60 

- 

- 

4 

- 

- 

- 

- 

- 

4 

10 

- 

5 

- 

- 

- 

- 

- 

- 

- 

1 

total 

1 

21 

124 

266 

266 

124 

21 

1 


Table 4. Minimal Betti table for the AGR example with 
n = 6, d = 5, and s = 24 


0 1 2 3 4 5 6 7 


0 

1 

- 

- 

- 

- 

- 

- 

- 

1 

- 

4 

- 

- 

- 

- 

- 

- 

2 

- 

32 

150 

256 

220 

96 

17 

- 

3 

- 

17 

96 

220 

256 

150 

32 

- 

4 

- 

- 

- 

- 

- 

- 

4 

- 

5 

- 

- 

- 

- 

- 

- 

- 

1 

total 

1 

53 

246 

476 

476 

246 

53 

1 


Table 5. Minimal Betti table for the AGR examples with 
n = 6, d = 5, and s > 42 


0 1 2 3 4 5 6 7 


0 

1 

- 

- 

- 

- 

- 

- 

- 

1 









2 

- 

56 

189 

216 

- 

- 

- 

- 

3 

- 

- 

- 

- 

216 

189 

56 

- 

4 









5 

- 

- 

- 

- 

- 

- 

- 

1 

total 

1 

56 

189 

216 

216 

189 

56 

1 


Table 6. The minimal Betti tables for the AGR examples 
with n = 6, d = 5, and s > 28 are the sum of 
Table [5] and the table below. 

0 1 2 3 4 5 6 7 

0 : - - - - - - - 

2: - - a(s ) fi(s) 7 (s) /3{s) a(s) - 

3: - a(s) /3(s) 7 (s) f3(s) a(s) 

4: - - - - - - - 

5: - - - - - - - - 

where a(s) := max{0, 6(31.5 — s)}, j3(s) := max{0, 15(36 — s)}, 
7 (s) := max{0, 20(42 — s)}. 
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Table 7. Non-minimal Betti table for the AGR examples 
with n = 6, d = 5, and s > 34 



0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

1 

- 

- 

- 

- 

- 

- 

- 

_ 

56 

210 

336 

280 

120 

21 

_ 

3 

- 

21 

126 

315 

420 

315 

126 

21 

4 

- 

6 

36 

90 

120 

90 

36 

6 

5 

- 

1 

6 

15 

20 

15 

6 

1 

total 

1 

84 

378 

756 

840 

540 

189 

28 


Table 8. Timings for the AGR examples with d = 5 and 
large s (in sec.) 


n 

Macaulay2 

Singular 

res 

Ires() 

LiftTree 

betti() 

minres () 

5 

2.47 

1.17 

0.06 

0.06 

0.49 

6 

87.74 

48.99 

0.73 

2.38 

23.39 

7 

1,782.44 

1,027.15 

9.97 

60.14 

322.45 

8 

50,203.20 

32,647.20 

137.25 

1,114.14 

9,002.31 

9 

- 

- 

1,852.42 

16,100.67 

- 

10 

- 

- 

24,281.08 

- 

- 


Table 9. Statistics for the non-minimal resolutions of the 
AGR examples with d = 5 and large s 


n 

sec. 

Terms 

#MuIt. 

#Add. 

#Canc. 

Qsparse 

5 

0.06 

59,903 

101,264 

44,790 

63 

0.282 

6 

0.73 

294,746 

447,232 

162,989 

327 

0.195 

7 

9.97 

1,292,567 

1,761,229 

496,922 

1,185 

0.130 

8 

137.25 

5,179,579 

6,433,983 

1,323,234 

3,562 

0.084 

9 

1,852.42 

19,311,659 

22,322,538 

3,166,754 

9,013 

0.053 

10 

24,281.08 

67,915,012 

74,543,926 

6,963,586 

20,518 

0.033 
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Figure 4. Grayscale image of the 756 x 840 matrix from 
the non-minimal resolution of the AGR example 
with n = 6, d = 5, and s = 42. Zero entries are 
white, entries with one term gray, and entries 
with two terms black. There are no entries with 



%/ /. 
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6.2. Canonical Nodal Curves and Prym Canonical Nodal Curves. 

Our next two series of examples are randomly constructed canonical nodal 
curves (CNC) (see Tables 10 and 0 and Prym canonical nodal curves 
(PCNC) (see Tables 12 and 13) of genus g. For the series of CNC examples, 


we work over the finite field with 32,003 elements, whereas we consider the 
PCNC examples over fields of random positive characteristic between 10,000 


and 30,000 as listed in the second column of Table 12 


Canonical curves are a widely studied topic, see, for example, [BS] for 
some references. Our interest in the minimal resolution of Prym canoni¬ 
cal curves comes from |CEFSj . where an analogue of Green’s conjecture for 
Prym curves was formulated. In terms of the resolutions computed with our 
method, both conjectures say that the degree zero parts of the syzygy matri¬ 
ces in the resolution have maximal rank. Computational results presented 
in |CEFSj indicate that the Prym-Green conjecture is very likely false for 
genus g = 8 and g = 16. It would be interesting to check experimentally 
whether g = 24 is an exception as well. But our timings suggest that this 
is way out of reach even for our improved method. In fact, computing a 
resolution with LiftTree for the case g = 24 would take, as a very rough 
estimate, about 100,000 years and 100 TB of memory. 

Remark 6.1. The non-minimal Betti numbers for the PCNC example with 
g = 14 as computed with LiftTree are shown in Table 14 Note that 


the fifth differential fa contains a 1932 x 1932 submatrix with constant 
entries. A black-and-white image of this submatrix is shown in Figure [5] 
Using specialized software such as FFLAS/FFPACK jDGPj which is eas¬ 
ily available via Sage [S + 14j . checking that this matrix has full rank takes 
only 0.92 seconds. Thus, including the time to compute the non-minimal 


resolution (224.15 seconds, cf. Table 12), the verification of the Prym-Green 


conjecture for this case takes 225.07 seconds using our new method. If we 
only consider the time to compute the non-minimal resolution up to the fifth 


differential (129.87 seconds, cf. Table 15), this can be reduced to 130.79 sec¬ 
onds. 

Remark 6.2. For the PCNC example with g = 16, the non-minimal Betti 
numbers as computed with LiftTree are shown in Table 16 Here, the 


sixth differential fa contains a 8910 x 8910 submatrix with constant entries. 
Using FFLAS/FFPACK [ DGP j. it takes 55.24 seconds to check that the 
kernel of this submatrix is one-dimensional. Including the time to compute 
the non-minimal resolution up to the sixth differential (10,776.54 seconds, 


cf. Table 18), it thus takes little more than 3 hours to check that the verifi¬ 


cation of the Prym-Green conjecture based on nodal curves fails in this case. 
This is a substantial improvement of the running time for this problem in 
comparison to the time needed in jCEFSj . 

We also managed to compute the minimal Betti numbers for this example 
using the Singular command bettiO, cf. Table 17 This computation 


took more than 37 hours, but it was run on a different compute server than 
the other examples due to the memory consumption exceeding 16 GB. 

Finally, starting from the non-minimal resolution, we even succeeded to 
compute the syzygy scheme of the Prym-Green extra syzygy which prevents 


the minimal resolution of the curve from being pure, cf. Table 17 However, 
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we found that in our example, this syzygy scheme is identical to the curve 
itself and does therefore not contribute any geometric information. 


Table 10. Timings for the CNC examples (in sec.) 


9 

Macaulay2 

Singular 

res 

Ires() 

LiftTree 

betti () 

minres () 

5 

0 

0 

0 

0 

0 

6 

0.01 

0 

0 

0 

0 

7 

0.07 

0.02 

0.01 

0 

0 

8 

1.41 

0.33 

0.02 

0 

0.04 

9 

25.45 

6.65 

0.07 

0.02 

0.97 

10 

452.10 

113.48 

0.31 

0.38 

9.28 

11 

7,945.04 

2,017.32 

1.83 

6.87 

155.52 

12 

- 

30,495.55 

15.58 

89.66 

1,238.28 

13 

- 

- 

142.20 

1,005.98 

23,877.02 

14 

- 

- 

1,351.59 

9,645.67 

- 

15 

- 

- 

12,935.45 

- 

- 


Table 11. Statistics for the non-minimal resolutions of the 
CNC examples 


9 

sec. 

^Terms 

#Mult. 

#Add. 

^Canc. 

Q sparse 

5 

0 

258 

466 

226 

0 

3.583 

6 

0 

1,411 

2,161 

831 

0 

2.520 

7 

0.01 

6,038 

8,104 

2,323 

3 

1.677 

8 

0.02 

22,343 

27,123 

5,453 

18 

1.075 

9 

0.07 

75,054 

84,804 

11,319 

63 

0.669 

10 

0.31 

235,179 

253,212 

21,434 

169 

0.408 

11 

1.83 

699,758 

730,490 

37,790 

380 

0.244 

12 

15.58 

1,998,583 

2,047,201 

62,944 

758 

0.144 

13 

142.20 

5,522,774 

5,593,998 

100,053 

1,393 

0.084 

14 

1,351.59 

14,854,811 

14,949,653 

152,971 

2,382 

0.049 

15 

12,935.45 

39,056,118 

39,164,376 

226,312 

3,866 

0.028 
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Table 12. Timings for the PCNC examples (in sec.) 


9 

Char. 

Macaulay2 

Singular 

res 

Ires() 

LiftTree 

betti() 

minres () 

6 

22,669 

0 

0 

0 

0 

0 

7 

10,151 

0.02 

0.01 

0 

0 

0 

8 

15,187 

0.18 

0.05 

0.01 

0 

0.01 

9 

18,947 

2.96 

0.75 

0.04 

0 

0.11 

10 

13,523 

64.78 

19.56 

0.15 

0.02 

1.96 

11 

25,219 

901.14 

344.62 

0.64 

0.40 

23.11 

12 

11,777 

16,597.90 

6,261.13 

3.43 

6.97 

305.45 

13 

24,379 

- 

- 

25.95 

98.34 

3,032.23 

14 

16,183 

- 

- 

224.15 

1,084.88 

40,106.95 

15 

20,873 

- 

- 

2,002.51 

10,953.23 

- 

16 

12,451 

- 

- 

18,612.82 

- 

- 


Table 13. Statistics for the non-minimal resolutions of the 
PCNC examples 


9 

sec. 

Terms 

#Mult. 

#Add. 

#Canc. 

Qsparse 

6 

0 

501 

534 

90 

0 

2.088 

7 

0 

2,818 

5,297 

2,663 

1 

2.271 

8 

0.01 

12,974 

24,165 

11,642 

8 

1.917 

9 

0.04 

48,711 

79,096 

31,409 

14 

1.390 

10 

0.15 

165,346 

232,455 

69,389 

43 

0.945 

11 

0.64 

524,473 

654,596 

135,098 

110 

0.617 

12 

3.43 

1,582,334 

1,812,443 

240,603 

285 

0.392 

13 

25.95 

4,594,249 

4,974,596 

401,889 

594 

0.244 

14 

224.15 

12,931,450 

13,525,642 

637,340 

1,154 

0.149 

15 

2,002.51 

35,482,705 

36,367,579 

970,090 

2,018 

0.090 

16 

18,612.82 

95,281,070 

96,541,345 

1,427,327 

3,409 

0.054 


Table 14. Non-minimal Betti table for the PCNC example 
with g = 14 



0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

0 : 

1 : 

1 

52 

303 

882 

1596 

1932 

1602 

903 

332 

72 

7 

“ 

2 : 

- 

17 

167 

738 

1932 

3318 

3906 

3192 

1788 

657 

143 

14 

total: 

1 

69 

470 

1620 

3528 

5250 

5508 

4095 

2120 

729 

150 

14 
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Table 15. Timings (in sec.) and statistics for the individual 
differentials fa in the non-minimal resolution of 
the PCNC example with g = 14. Note that fa 
is given by the input ideal. 


i 

^Generators 

fa Terms 

Qsparse 

Time 

~ fa Terms/sec. 

1 

69 

3,064 

44.406 

- 

- 

2 

470 

50,208 

1.548 

0.54 

92,978 

3 

1,620 

371,814 

0.488 

6.25 

59,490 

4 

3,528 

1,290,516 

0.226 

37.98 

33,979 

5 

5,250 

2,639,132 

0.142 

85.10 

31,012 

6 

5,508 

3,436,908 

0.119 

70.07 

49,050 

7 

4,095 

2,917,668 

0.129 

21.21 

137,561 

8 

2,120 

1,593,830 

0.184 

2.72 

585,967 

9 

729 

530,548 

0.343 

0.25 

2,122,192 

10 

150 

94,488 

0.864 

0.02 

4,724,400 

11 

14 

6,338 

3.018 

0.01 

(633,800) 

2 to 11 

23,484 

12,931,450 

0.149 

224.15 

57.691 


Table 16. Non-minimal Betti table for the PCNC example 
with g = 16 



0 

1 

2 

3 

4 

5 

6 

7 

8 

0 : 

1 

- 

- 

- 

- 

- 

- 

- 

- 

1 : 

- 

75 

539 

1980 

4653 

7590 

8910 

7623 

4730 

2 : 

- 

19 

225 

1221 

4015 

8910 

14058 

16170 

13662 

total: 

1 

94 

764 

3201 

8668 

16500 

22968 

23793 

18392 





9 

10 

11 

12 13 







2079 

615 

110 

9 - 







8415 

3685 

1089 

195 16 







10494 

4300 

1199 

204 16 



Table 17. 

Minimal Betti table for the PCNC example with 




9 = 16 







0 

1 

2 

3 

4 

5 

6 

7 

8 

0 : 

1 

- 

- 

- 

- 

- 

- 

- 

- 

1 : 

- 

75 

520 

1755 

3432 

3575 

1 

- 

- 

2 : 

- 

- 

- 

- 

- 

1 

6435 11440 

11583 

total: 

1 

75 

520 

1755 

3432 

3576 

6436 11440 

11583 





9 

10 

11 

12 13 




7800 3575 1080 195 16 


7800 3575 1080 195 16 
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Table 18. Timings (in sec.) and statistics for the individual 
differentials fa in the non-minimal resolution of 
the PCNC example with g = 16. Note that fa 
is given by the input ideal. 


i 

^Generators 

$Terms 

Qsparse 

Time 

~ fa Terms/sec. 

1 

94 

4,706 

50.064 

- 

- 

2 

764 

96,185 

1.339 

1.75 

54,963 

3 

3,201 

881,460 

0.360 

51.64 

17,069 

4 

8,668 

3,884,192 

0.140 

669.57 

5,801 

5 

16,500 

10,428,510 

0.073 

3,332.74 

3,129 

6 

22,968 

18,632,465 

0.049 

6,720.84 

2,772 

7 

23,793 

23,035,835 

0.042 

5,663.46 

4,067 

8 

18,392 

19,959,667 

0.046 

1,897.21 

10,521 

9 

10,494 

12,027,377 

0.062 

260.48 

46,174 

10 

4,300 

4,889,758 

0.108 

14.46 

338,158 

11 

1,199 

1,257,155 

0.244 

0.59 

2,130,771 

12 

204 

178,565 

0.730 

0.08 

2,232,063 

13 

16 

9,901 

3.033 

0 

- 

2 to 13 

110,499 

95,281,070 

0.054 

18,612.82 

5,119 
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Figure 5. Black-and-white image of the 1932 x 1932 con¬ 
stant submatrix from the non-minimal resolution 
of the PCNC example with g = 14. Zero entries 
are white, constant entries are black. Note the 
somehow fractal structure. 
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